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ABSTRACT 

We report the discovery of the merging cluster, RXCJ2359.3-6042, from the REFLEX II cluster survey and present our results from all 
three detectors combined in the imaging and spectral analysis of the XMM-Newton data. Also known as Abell 4067, this is a unique 
system, where a compact bullet penetrates an extended, low density cluster at redshift z=0.0992 clearly seen from our follow-up 
XMM-Newton observation. The bullet goes right through the central region of the cluster without being disrupted and we can clearly 
see how the bullet component is stripped of its layers outside the core. There is an indication of a shock heated region in the east of 
the cluster with a higher temperature. The bulk temperature of the cluster is about 3.12(±0.13)keV implying a lower mass system. 
Spearheading the bullet is a cool core centred by a massive early type galaxy. The temperatures and metallicities of a few regions 
in the cluster derived from the spectral analysis supports our conjecture based on the surface brightness image that a much colder 
compact component at 1.55(±0.10)keV with large metallicity (0.75 Z Q ) penetrates the main cluster, where the core of the infalling 
component survived the merger and left stripped gas behind at the centre of the main cluster. We also give an estimate of the total 
mass within r 50 o, which is about 2xl0 14 M Q from the deprojected spherical /3 modelling of the cluster in good agreement with other 
mass estimates from the M-Tx and M-<x v relations. 
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1. Introduction 


Merging galaxy clusters have been found to be interesting 
laborat ories for the study of a variety of astrophysica l pro¬ 
cesses dFeretti et al .1120021 iMarkevitch & Vikhlininll2007l) . They 
are also the major process in which galaxy clusters grow and 
have therefore a large influence on the global properties of these 
systems. Observationally we see a wide range of different mor¬ 
phologies in cluster mergers, which is due to different merger 
geometries and the variety of morphologies of the merging clus¬ 
ters including different mass ratios. 

Based on the comple te flux-limited cluste r survey REFLEX 
II in the southern sky (Bohringer et alJl2013l) comprising more 
than 900 galaxy clusters, we constructed a volume-limited sub¬ 
sample (Bohringer et al., in prep.), in which we study 108 of 157 
clusters in XMM-Newton follow-up observations for a syst em¬ 
atic study of galaxy cluster str ucture statistics (Bohringer et alJ 
2010; Chon & Bohringerl l2012l) . 


In this survey we found a very interesting merging cluster, 
RXCJ2359.3-6042, also known as Abell 4067, at a redshift of 
0.0992 with a special merger configuration that has never been 
observed in this form before. In this system we witness the en¬ 
counter of two very dissimilar systems, a smaller compact sys¬ 
tem with a cool component, which penetrates a larger, more dif¬ 
fuse cluster with low central density in a bullet-like fashion. The 
fact that a compact bullet merges with an extended cluster with¬ 
out a dense core r esembles the case of IE 0657-558, known as 
the Bullet Cluster dMarkevitch et alJ d2004l) ). Unlike this cluster, 
however, the relative velocities and the Mach number associated 
with RXCJ2359.3-6042 are much lower, and so the interaction 
processes are quite different. In particular, the interaction seems 


to be less violent and probably as a result we clearly see a trail 
of debris behind the bullet which is not observed in the Bullet 
Cluster. This can provide interesting information on the strip¬ 
ping process acting on the infalling component. The mass of the 
cluster and the associated ICM temperature are also lower than 
in the Bullet Cluster making temperature measurements easier. 
These circumstances make this object unique as a laboratory for 
the study of the stripping process as explained below. 

Throughout the paper we adopt a flat A-cosmological model 
with Q m =0.3 and a Hubble constant of Ho=70 kms^'Mpc -1 for 
the calculations of all distant dependent quantities. We used Nh 
of 1.41 x 10 2,) cirr 2 for the spectral analysis. All quoted uncer¬ 
tainties refer to one sigma error. 


2. XMM-Newton observations 

The digitised sky-survey (DSS) optical image is shown with the 
XMM-Newton X-ray contours in Fig. Q] which gives a first im¬ 
pression of the geometric configuration of RXCJ2359.3-6042. 
On top of the extended X-ray emission from the cluster we 
see a smaller, much stronger emission from a source, which 
coincides with an optical galaxy, near the centre of the image 
with an extended tail towards the south-west. This galaxy is at 
z=0.100 dTeague et al .11 1 9901) . 

We used the XMM-Newton observation of RXCJ2359.3- 
6042 with Obs.ID 0677180601 as a part of the REFLEX II clus¬ 
ter survey project to investigate further. After the flare cleaning 
both MOS have 11 ks exposure while pn is left with 8 ks. Stan¬ 
dard data analysis was performed with SAS vl3.0.1. Our basic 
X-ray i maging data reduction follows closely the procedures out¬ 
lined in iChon et al. ( 2012 ). and we briefly summarise them here. 
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the values of all background components were fixed to the pre¬ 
determined values from the background region except for the 
normalisations, which are scaled appropriately according to the 
areas. There are known cross-calibration issues among the EPIC 
cameras (iRead et alJ 2014); h owever, this is less of a problem for 
cooler clusters (see e.g. ISchellenberger et al.1 (12014]) ). Hence, in 
the final fitting the temperature, metallicity, and normalisation 
of the cluster were kept free, and all three detectors were linked 
apart from the normalisation. 


3. Interpretation of the data 

3.1. Image analysis 

Figure [2] shows a combined flux image of three XMM-Newton 
detectors of the cluster in the [0.5-2] keV band, smoothed with 
4” kernel. 


Fig. 1. RXCJ2359.3-6042 shown in the DSS optical image overlaid 
with the XMM-Newton X-ray contours in the [0.5-2] keV band from 
the same data as in Fig. [2] The contours are obtained from a background 
subtracted image where the first contour traces 3x 10~ 5 counts s _1 pix~‘. 
The flux of the contours increases by steps of 1.5 subsequently. The 
compact component with high X-ray emission is centred at the galaxy 
at 23 ,! 59"'04.2’ in R.A. and -60‘ , 36"'34 I in Dec. 


The XMM-Newton observation for all three detectors were flare- 
cleaned and out-of-time foot) events were statistically subtracted 
from the pn data. Point sources were removed, and the modelled 
background was subtracted b a sed on the blank sky observation 
provided bv iRead & Ponmanl ( 2003 ). The cleaned images from 
all three detectors were combined with an exposure correction. 

For the spectral analysis we used Filter Wheel Closed (FWC) 
data to remove the contribution from the particle background. 
The FWC spectrum was scaled according to the ratio of the cor¬ 
ner events between the observation and FWC data. We also re¬ 
moved the oot events of the FWC-pn data before scaling, which 
mostly affects the very low energy rangcQ The region outside 
rsodH °f the cluster within the field of view was taken as a 
background from which the contributions from the X-ray cos¬ 
mic background together with the residual particle background 
were determined. We considered an unabsorbed thermal model 
for the Local Hot Bubble, unresolved point sources described 
by an absorbed power law, and a cool absorbed thermal halo. 
The unabsorbed thermal component is best fitted with 0.09 keV 
while the absorbed thermal model with 0.29 keV. In both cases 
the metallicity is set to unity and the redshift to zero. The power 
law describing unresolved point sources is best fitted with an in¬ 
dex of 1.41. The residual particle background was modelled with 
a power law, where the best fit indices for MOS1, MOS2 and 
pn are 0.1, 0.41, and 0.46. The spectral regions were then fitted 
with a thermal component representing the cluster emission to¬ 
gether with the background components described above using 
the APEC model in XSPEC. Given the limited photon statistics 



Fig. 2. XMM-Newton flux image of RXCJ2359.3-6042. All three de¬ 
tectors were combined after the exposure correction in the [0.5-2.0] keV 
band. A scale of three arcmin is shown with a bar. 


The main cluster shows a very shallow, extended X-ray sur¬ 
face brightness distribution. It appears as a dynamically young 
system which has not virialised, and has not formed a dense 
core. Alternatively, one can attribute the shallow surface bright¬ 
ness distribution to the disturbance caused by the merger itself, 
but in any case the main cluster could not have been a very regu¬ 
lar cluster with a dense core before the merger, as we will argue 
below. Embedded in this low surface brightness cluster we see 
a small, bright compact component with a trail of patchy X-ray 
emission behind. We extracted a radial surface brightness profile 
of this component, shown in Fig.[3]as a solid line. This profile is 
calculated only from the leading edge, which is almost identical 
to that of the entire region of the compact component. This con¬ 
firms that the extended emission does not result from the debris 
trailing from the compact component. The centre of the profile 
is taken at the maximum of the X-ray emission, which coincides 


‘Newly added FWC oot pn data are found at 
http: //xmm2 . esac. esa. int/external/xmm_sw_cal/background/f'fl'ifeiL!6iyRtfH] , [)ii7 i; ?'t t ^n , ^F? l i_ I 2 , 0 l f'4_4 1 l PSfrtmF- Clearly it shows 
' The radius within which the mean cluster mass density is 500 times extended emission, which can be compared to the radial pro- 
the critical density of the Universe. To calculate^ for RXCJ2359.3- hies of typical point sources in the same observation shown as 

6042 we used the recipes as described in lBohringer et ali ( 20131 1. dashed and dash-dotted lines. The most obvious interpretation of 
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distance [arcmin] 

Fig. 3. Profiles of the cool and extended component penetrating the 
cluster (solid) in the [0.5-2.0] keV band in comparison to two point 
sources in the same observation. One point source is located at a similar 
off-axis angle (dash-dotted) and the other point source is at about twice 
the off-axis angle (dashed) as the cool component. 

Fig. [2]is then that a compact cool core cluster has fallen into the 
diffuse, extended main system, and the trail of gas seen in the 
middle of the picture is a mixture of the gas stripped from the 
infalling bullet and main cluster. Inspection of the optical image 
in Fig.Q]shows that the bright dominant galaxy is still at the cen¬ 
tre of the infalling compact cluster. With the 18 galaxy redshifts 
available in the NEE@ within the rgoo of RXCJ2359.3-6042 we 
infer the approximate direction of the infall. For a cluster mass 
of 2xlO 14 M 0 , as will be estimated in Sect. 4, the expected infall 
velocity is about 1310 k m s 1 . The line-of-sight relative velocity 
difference between the infalling bright core with respect to the 
cluster redshift is about 430 km s~ 1 . This implies that the infall 
direction is close to the plane of the sky within about 20 degrees. 
This configuration allows a relatively simple interpretation of the 
merging process limiting the ambiguity of projection effects. 

The aforementioned features reveal a particular case that the 
“bullet” has gone through the centre of the cluster, and still has 
preserved its core. In this case we expect the metallicity of this 
core region to be large, which indica tes that it was undisturbed 
for long time (iBohringer et al.ll2004l) . From this we can further 
infer that main cluster did not have a dense core itself, otherwise 
the collision of the two cores would have disrupted the infalling 
core. We can further deduce that if the main cluster did not have 
a cool core, the patchy emission that we observe as a trail of the 
“bullet” must be the debris of the outer layers of the infalling 
system. 

Further inspection of the image shows a sharper edge of the 
cluster on the east side than on the western counterpart, and faint 
emission outside the east edge. Analogous to the bullet cluster, 
we expect the sharp edge (approximately bounded by the west¬ 
ern sides of the green rectangular regions in Fig. [4} to be the 
contact discontinuity of the bullet material. In this case we would 
expect to see a shock further to the east of this edge, but because 
of the low surface brightness and poor photon statistics we can¬ 
not detect a surface brightness jump that would reveal the loca¬ 
tion of this shock. Nevertheless, as another consequence of such 
a shock, we would expect the shock heated region to the east of 

:,: The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, un¬ 
der contract with the National Aeronautics and Space Administration. 



Fig. 4. Defined regions for the spectral analysis. The central 0.7xr 50 o 
is marked by a green circle centred at the yellow cross along with 
three boxes that divide the main cluster. The centre of the cluster is 
at 23 a 58'"52.3 ! in R.A. and -60 rf 37" ! 59 I in Dec. The cool component 
(CC) is marked by a red circle centred at the highest X-ray emission. 
The suspected shocked regions are enclosed by two green boxes. 


the contact discontinuity to be quite hot. An obvious difference 
to the Bullet Cluster is that we have a very shallow opening an¬ 
gle of the cone of the contact discontinuity which is expected to 
follow the shape of the possible Mach cone. This implies a very 
mild shock. The possible shock wave still needs to be located 
and needs to be better quantified together with the cold front 
with better data in the future. 

3.2. Spectral analysis 

In this section we will study the temperature and metallicity dis¬ 
tribution in the cluster, which will help to verify our speculations 
about the history of the merger laid out in the previous section. 

Figure |4] defines the regions from which we extracted the 
spectra. The cluster bulk properties are determined from the X- 
ray emission inside the large circle enclosing 0.7xr5oo excluding 
the circular region centred at the infalling core, which was anal¬ 
ysed separately. Given the extent of the cool component (CC) 
found in the surface brightness profile, the CC region was de¬ 
fined by a radius of one arcmin. The main part of the cluster 
is divided into three boxes; Boxes 1 and 3 are considered to be 
less affected by the infalling component than Box 2. We also ex¬ 
cluded the CC from Box 2 when extracting the spectrum. The 
suspected shock-heated regions are marked by two green rect¬ 
angles. The corresponding best fit temperatures and metal abun¬ 
dances for these regions are given in Table Q] 

The temperature measured within 0.7xr5oo is 3.12keV with 
a metallicity of 0.17Z o without the CC. The bright, compact 
CC is found to be at the lowest temperature in the system with 
the largest metal abundance. This nicely confirms our conjecture 
about this region being an old cool core. The next coldest com¬ 
ponent is found in Box 2, which we assumed to be the debris of 
the infalling system mixed with the main cluster. Boxes 1 and 
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region 

T (keV) 

z Q 

CC 

1.55 + 0.10 

0.75 + 0.19 

cluster (-CC) 

3.12 + 0.13 

0.17 + 0.04 

Boxl 

3.54 ± 0.26 

0.18 + 0.09 

Box2 (-CC) 

3.08 + 0.12 

0.27 ± 0.08 

Box3 

3.51 ±0.22 

0.18 + 0.07 

shock (N+S) 

4.5 ± 0.59 

0.35 ± 0.22 


Table 1. Best fit temperature and metallicity of the spectral areas de¬ 
fined in Fig.[4]with lcr errors. 


3 contain relatively hot gas in comparison. In particular, we find 
the region marked in green rectangles on the eastern front, which 
we assume to be shock-heated, to be the hottest region in the 
cluster. This is analogous to the shock heated region in the bul¬ 
let cluster, but everything is scaled to lower temperature. Hence, 
the temperature and metallicity distributions support our conjec¬ 
ture about the merging process laid out in the previous section. 
Currently, we are limited by the photon statistics and can neither 
pin-point the exact location of the shocked region nor the contact 
discontinuity with the spectral analysis. 

4. Modelling the main cluster component 

To obtain more insight into the properties of the main cluster 
component, we use a spherical approximation to the shape of the 
cluster, excising a circular area around the CC and point sources. 
We fit a /j-model with a constant background to the surface 
brightness distribution of the cluster whose centre, marked by 
a yellow cross in Fig. [4] is determined by the centre shift method 
described in lChon et alJ (120121) . The best fit value of /j is found 
to be 0.64+0.062, and that of the core radius to be 3.04+0.32 
arcmin, both with lcr error. A variation of ~10% in the back¬ 
ground corresponds to a change in the estimates of two parame¬ 
ters similar to the errors. With these values and the cluster bulk 
temperature found in the previous section we can estimate the to¬ 
tal mass of the main cluster assuming a spherical symmetry and 
isothermality of the ICM, which allows an analytical deprojec¬ 
tion of the observed surface brightness. We then get an estimate 
of the total mass of 2.1(±0.37)xl0 14 M o for the cluster temper¬ 
ature of 3.54keV. We consider the temperatures in Boxes 1 and 
3 of 3.54keV a better representation of the bulk temperature of 
the undisturbed main cluster, compared to the temperature mea¬ 
surement in Box 2 or the total cluster region, since the mixture 
of main cluster gas with colder gas stripped away from the bullet 
will bias the single temperature fit low. Nevertheless, to bracket 
the uncertainties in the temperature estimate, we also estimated 
the mass of the cluster for two temperature limits, 3 and 3.8 keV. 
For this temperature range we find that the estimated mass lies 
between 1.8 and 2.2xl0 14 M o . 

It is also interesting to compare this mass estimate to other 
mass probes. With the estimated cluster temperature we can use 
the M-Tx scaling relation, for example bv lArnaud et alJ d2005l) . 
to obtain the cluster mass. Using the overdensity 6 of 500 we get 
a mass estimate of 2(±0.25)xl0 14 M o for the cluster temperature 
of 3.54keV. For 3.12keV we get a roughly 5.5% lower mass, 
hence in both cases we have a good agreement with our hydro¬ 
static mass estimate. We will consider another probe of mass in 
the next section with optical spectroscopic data. 

From the /j-model fit to the cool core excised main 
cluster image, we determine a central electron density of 
about 2.5xl0~ 3 cm~ 3 . This yields a central entropy value 
of 190+14 keV cm' 2 , which is well in the regime of non¬ 


cool core clusters (see iPratt et all ( 2010 11. It is less ex¬ 
treme than that o f Coma with its central entropy value 
of 535 keV cm~ 2 (iFusco-Femiano et alJ 1201 3l) . and that of 
the low surface brightness cluster A76 with a value of 
400 keV cm -2 (iOta et alil2013l ). but it shows that the properties 
of the main cluster is far from the stage where it could have de¬ 
veloped a cool core, in agreement with our conclusion above. 

4.1. Dynamical mass estimate 

An estimate of the dynamical mass of a cluster can be inferred 
from measurements of spectroscopic redshifts of member galax¬ 
ies. In the presence of a large number of redshifts of a hundred or 
more a detailed analysis of the redshift distribution is appropriate 
as outlined bv lBeers et ali d 19901) . for example. For smaller data 
sets with at least ten redshifts a scaling relation of velocity dis¬ 
persion and cluster mass as provided for example by Biviano et 
al. (2006) is considered a reasonable and robust approach. This 
calibration of the scaling relation is based on the study of 64 
galaxy clusters in cosmological N-body simulations. They con¬ 
sider that mass estimates are reliable for a sample that involves 
at least ten galaxy redshifts within the virial radius. There are 
19 spectroscopic redshifts available in the NED database for 
RXCJ2359.3-6042 within r^ pn, with w hich we calculated the ve¬ 
locity dispersion based on iHarrisonl dl974l) and iBiviano et alJ 
d2006l) . We applied a three-sigma clipping, which removes one 
galaxy, to obtain the final dispersion. We followed the recipes 
given in iBiviano et al.1 ( 2006 ) to calculate a mass estimate con¬ 
sidering their Eq. 2 and Fig. 4, and obtain a mass estimate 
of 2.2 x10 14 M o for the one-dimensional velocity dispersion of 
605+87 km s^ 1 and h-0.1. The quoted error on the velocity dis- 
person results from 1000 bootstrap realisations, which were sub¬ 
ject to the same procedure as were the data, which results in a 
43% error for the mass estimate including the uncertainty of the 
scaling relation given in IBiviano et al.l (i2006l) . Since the dynam¬ 
ical mass estimate is based on the dispersion of radial velocities, 
fair agreements between different mass estimates of this system 
implies that the dynamical disturbance cannot be large in the ra¬ 
dial direction. It provides another signature that the merger is 
close to the plane of the sky, in agreement with the angular tra¬ 
jectory of the bullet component that we calculated earlier. 


5. Summary and outlook 

We presented the analysis of the XMM-Newton data of a merg¬ 
ing system, RXCJ2359.3-6042, which shows a main cluster with 
a low X-ray surface brightness being penetrated by a cooler com¬ 
pact component whose core is still intact. Our conjecture about 
the merging configuration is supported by the spectral tempera¬ 
tures and metal abundances found in different regions of the clus¬ 
ter. We find that the mass of the main cluster without the infalling 
component is about 2xlO 14 M 0 , which is based on three differ¬ 
ent mass estimates for a bulk temperature of 3.5 keV. Within the 
current limit of photon statistics we can tentatively infer that the 
infalling component must have induced a mild shock given the 
wide opening angle. A deeper observation will offer a chance to 
quantify the processes ongoing in this special object in greater 
detail. In particular, it may allow us to identify the exact location 
of the contact discontinuity, to obtain more quantitative details 
about the shock, and to reconstruct some of the properties of the 
infalling system such as a limit to the total gas mass. 
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